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Abstract 

We provide an exact calculation of next-to-next-to-leading order (NNLO) massive corrections to Bhabha scattering 
in QED, relevant for precision luminosity monitoring at meson factories. Using realistic reference event selections, 
exact numerical results for leptonic and hadronic corrections are given and compared with the corresponding ap- 
proximate predictions of the event generator BabaYaga@NLO. It is shown that the NNLO massive corrections are 
necessary for luminosity measurements with per mille precision. At the same time they are found to be well ac- 
counted for in the generator at an accuracy level below the one per mille. An update of the total theoretical precision 
of BabaYaga@NLO is presented and possible directions for a further error reduction are sketched. 
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1. INTRODUCTION 

Because of its experimental and theoretical charac- 
teristics, Bhabha scattering is the prime process used at 
e + e~ colliders to monitor their luminosity |1|. At the 
GeV-scale e + e~ colliders, from <p to B factories, large 
angle Bhabha scattering is measured with experimen- 
tal uncertainties at the per mille level. Correspondingly, 
the Monte Carlo (MC) generators used in the data anal- 
ysis, like BabaYaga@NLO EE II, BHWIDE and 
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MCGPJ |6|, are precision tools that include the large 
logarithmic contributions due to soft and collinear mul- 
tiple photon emission matched with the exact next-to- 
leading order (NLO) photonic corrections. In addition 
to the these ingredients, the MC codes include the effect 
of leptonic and hadronic vacuum polarisation, the latter 
computed using a data based routine for the evaluation 
of the non-perturbative light quark contribution 171 . 

The above framework implies that the MC pro- 
grams are affected by a theoretical uncertainty stem- 
ming mainly from missing (subleading) contributions 
at the level of NNLO radiative corrections. Therefore 
exact NNLO calculations are the benchmark for a reli- 
able assessment of the error associated to the MC pre- 
dictions. Since the whole class of the NNLO QED cor- 
rections to Bhabha scattering became recently available 
(see 1 1 1 for a review), various comparisons between the 
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exact calculations and the MC results allowed to con- 
clude that the theoretical accuracy of the most precise 
luminosity tools is at the level of one per mille. Al- 
though this uncertainty seems sufficient by comparison 
with the present experimental error, there is still room 
for a more reliable assessment and, presumably, a fur- 
ther error reduction. 

Here we focus on the NNLO massive corrections to 
Bhabha scattering, for which only preliminary results 
limited to leptonic corrections for a few experimental 
selections were used in the most recent, official estimate 
of the MC precision [ 1 1. These limitations are overcome 
in the present study [8 1, wherein the full set of NNLO 
leptonic and hadronic corrections is computed and all 
the event selections of experimental interest at meson 
factories are addressed. The exact NNLO results are 
compared with the corresponding 0(a 2 ) truncated pre- 
dictions of the MC BabaYaga@NLO. 

2. THE NNLO LEPTONIC AND HADRONIC 
CORRECTIONS 

The full treatment of Nf = 1,2 massive corrections 
to Bhabha scattering requires the calculation of i) pure 
two-loop diagrams, including irreducible vertex and box 
corrections that are known to give rise to large lead- 
ing L?,L - ln(s/m 2 ), collinear contributions ii) loop- 
by-loop corrections Hi) real photon corrections and iv) 
corrections due to real lepton and hadron pair emis- 
sion. The latter generate L? contributions canceling 
those from irreducible virtual loops. 

2.1. Exact calculation 

Exact electron loop corrections (Nf = 1) were first 
computed in ll9l [T0ll . while heavy fermion and hadronic 
loop contributions (Nf = 2) were later calculated in 
lUIlEIinilllllTSl. As said above, these virtual con- 
tributions must be combined with the real corrections. 
The cross section with complete NNLO leptonic and 
hadronic corrections can be written as 

J NNLO J^-NNLO Hrr NLO (f,{\ 

aCr N f =l,2 _ d<T vin "(TysoU^) 

dQ ~ dQ + dQ 

dcr NL ° .(ai) dcr L ° 

y.hai-d v ' uu re al 

dQ dQ. K ' 

which can be rearranged as 

dQ dQ + dQ + dQ 

where w is a soft-hard photon separator. The terms 
entering Eq. |2ji have been computed using various 



semi-analytical and MC tools developed over the years, 
i.e. by the package bha_nnloJif |[T6l for the vir- 
tual + soft corrections dcr™^°(to)/a?Q, the MC code 
BHAGEN-1PH Ell QD for hard photon radiation 
dcr^^^/dQ, and the generators HELAC-PHEGAS 
IITgi ' and EKHARA (20] El for real lepton and hadron 
pair emission dcr^/dQ, respectively. For the R 
contribution to the hadronic vacuum polarisation we 
adopt the recent compilation encoded in the routine 
VP_HLMNT_v2_0 (H |23] |24]]. Strictly speaking, 
EKHARA allows the computation of real pion pair 
emission only and no other generator exists for the cal- 
culation of similar processes containing other hadron 
pairs. However, we observed in our numerical study |8] 
that the cross section for real pion pair emission is al- 
ways very small at all meson factories (at the level of 
10~ 5 times the Born cross section or below it). There- 
fore an educated guess is that the processes of pair emis- 
sion containing hadrons heavier than pions have a neg- 
ligible effect in our numerical predictions. 

2.2. Exact numerical results 

In Table Q] we show the relative contribution of the 
exact NNLO massive corrections for all the interesting 
experiments at meson factories and using realistic refer- 
ence event selections used in luminosity measurements. 
In Table[I]S A = o-f LO /o- BY with x = e + e', lep, had, tot, 
where tot stands for the sum of leptonic (lep) and 
hadronic (had) corrections, and <x BY is the full cross sec- 
tion from BabaYaga@NLO. We refer to |8| for details 
about the cuts used in our analysis. 

Some comments are in order here. 

Table 1: Exact relative NNLO massive corrections S x = o"? NLO /o"BY: 
in per mille. 

S lep Shad S tot 

KLOE 1.020 -3.935(5) -4.472(5) 1.02(4) -3.45(4) 

BES 3.097 -2.246(8) -2.771(8) 

BES 3.650 -1.469(9) -1.913(9) -1.3(1) -3.2(1) 

BES 3.686 -1.435(8) -1.873(8) 

BaBar 10.56 -1.48(2) -2.17(2) -1.69(8) -3.86(8) 

Belle 10.58 -4.93(2) -6.84(2) -4.1(1) -10.9(1) 



According to the above discussion, the results for the 
hadronic corrections have been obtained using the full R 
parameterization for the virtual corrections and the con- 
tribution of pions only for real pair emission. Moreover, 
we do not provide results for the hadronic corrections 
at the BES (center of mass) cm. energies on top of the 
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Table 2: Comparison of the exact massive NNLO with BabaYaga@NLO results, in per mille. 
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J /if/ and i//{2S) resonances. Actually, we observed JU 
that for such energies the NNLO Bhabha cross section 
is dominated by the contribution of narrow resonances 
that can not be treated like mere perturbative effects and 
require a separate study. 

Among the corrections induced by the different par- 
ticle species, it can be seen from Table [T] that the domi- 
nant contribution is given by the electron pairs, with an 
increasing importance of muons and hadrons at the B 
factories, the t contribution being always negligible at 
meson factories. 

In particular, among the real emission processes we 
concluded that only the reaction <? + e~ — > e + e~e + <?~ gives 
significant contributions to the cross section used in the 
luminosity measurements. When the accuracy of the ex- 
periment reaches the per mille level, this process has to 
be considered and its contributions added to the theoret- 
ical cross section or alternatively subtracted as a back- 
ground from the experimental cross section. 

The main conclusion of the exact NNLO calculation 
is that the total correction is about 0.3 - 0.4% at KLOE, 
BES and BaBar and reaches about 1% at Belle. There- 
fore, complete NNLO massive corrections, or a their ap- 
proximation, are certainly needed for precise luminosity 
calculations. 

2.3. The NNLO massive corrections in the code 
BabaYaga@NLO 

BabaYaga@NLO is one of the most precise and 
widely used theoretical tools to monitor the luminosity 
of meson factories. It is a QED Parton Shower gen- 
erator, whose intrinsic leading logarithmic accuracy is 
improved by the inclusion of exact NLO soft+virtual 
and hard photon correction factors. In BabaYaga@NLO 
these NLO factors are dressed by self-energy insertions, 



so that a subset of the complete NNLO massive cor- 
rections previosuly discussed is included in the code. 
More precisely, the contributions accounted for are: i) 
the factorizable loop-by-loop corrections within the full 
class of NNLO virtual corrections ii) the contribution of 
real (soft and hard) photon emission in the sector of real 
corrections. In other words, in the generator any con- 
tribution from real pair radiation and irreducible virtual 
corrections is neglected, which makes the approach the- 
oretically consistent because it avoids an imbalance of 
L? contributions. 

Both leptonic and hadronic self-energy contributions 
are taken into account in the code. For the R pa- 
rameterization we use in the present study the same 
VP _HLMNT_v2_0 routine as in the exact calculation, in 
order to perform consistent comparisons. 

3. EXACT NNLO MASSIVE CORRECTIONS VS. 
BABAYAGA @ NLO 

The quality of the approximation inherent in 
BabaYaga@NLO is shown in Table [2] in comparison 
with the results of the exact calculation. Again realis- 
tic reference event selection cuts are considered. 

It can be seen from Table [2] that there is a gener- 
ally good agreement for all the meson factories exper- 
iments, the relative differences always being below the 
one per mille. There is a maximum registered difference 
of about 0.07% at Belle, where the corrections have the 
largest impact, but the agreement is at the level of a few 
units in 10 4 for all the other experiments. Therefore the 
very bulk of NNLO massive corrections is accounted for 
in BabaYaga@NLO. 

In our analysis, we also studied how the above con- 
clusions are stable against variation of the experimental 
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Figure 1: The relative difference, in per mille, of the NNLO mas- 
sive leptonic and hadronic corrections between exact cr^^f and 
BABAYAGA@NLO cr N B ^ L0 ', as a function of an acollinearity cut for 
two different angular acceptance regions for KLOE-like event selec- 
tions. 
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cuts, such as acceptance and acollinearity cuts. We no- 
ticed that for the leptonic corrections there is a slight in- 
crease in the difference for particularly tight acollinear- 
ity cuts but in general the difference is not particu- 
larly sensitive to the variation of the event selection. 
In any case, also for particularly severe cuts, the max- 
imum observed difference still remains below the one 
per mille [8]. An example of such a cut-dependence 
study is shown in Fig. [T]for KLOE, which is the exper- 
iment with the smallest uncertainty in the measurement 
of the luminosity, at the level of 0.2%. It can be seen 
that a variation of the acollinearity cut from its refer- 
ence value of 9° to smaller values only slightly affects 
the relative difference between the exact calculation and 
the BabaYaga@NLO approximation for leptonic cor- 
rection, leaving the situation for hadrons practically un- 
changed. 

4. CONCLUSIONS 

We provided an exact calculation of the NNLO mas- 
sive corrections to Bhabha scattering in QED, exploit- 
ing a number of analytical results and computational 
tools developed over the years. We presented exact 
numerical predictions for leptonic and hadronic cor- 
rections in the presence of realistic cuts for luminos- 
ity monitoring at meson factories. We compared the 
exact results with the approximate ones provided by 
the MC code BabaYaga@NLO. We concluded that 
NNLO massive corrections are well accounted for in 
the generator at an accuracy level below the one per 



mille. This reinforces the estimate of the total theo- 
retical precision of BabaYaga@NLO, previously estab- 
lished on the grounds of only partial results for NNLO 
massive corrections. The present total accuracy of 
BabaYaga@NLO can be estimated 0.1% at KLOE and 
BES (excluding data taking on top of narrow resonances 
J/if/ and ip(2S)), and about 0.15% at the B factories. 

This precision is today sufficient by comparison with 
the experimental accuracy of meson factories. Nonethe- 
less, there is still room for a further error reduction 
along the following directions. There is a need for 
an assessment of the theoretical accuracy for luminos- 
ity measurements in a close vicinity of narrow reso- 
nances. This requires detailed studies of the uncertain- 
ties associated to the hadronic vacuum polarisation and 
NNLO hadronic corrections, including beam spread ef- 
fects. The leading part of the presently missing NNLO 
massive corrections in BabaYaga@NLO could be im- 
plemented through e.g. QED Structure Functions, to ac- 
count for the interplay between virtual irreducible cor- 
rections and lepton pair radiation. Last but not least, the 
recently obtained exact results for the one-loop correc- 
tions to radiative Bhabha scattering l25l . as well as for 
<? + e~ — > y^iCy [25 26|, are the final NNLO benchmark 
that should be deeper investigated in comparison with 
the present MC approximation. 

All these perspectives are left to future works. 
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